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A spin-glass transition occurs both in and out of the limit of validity of mean-field theory on 
a diluted one dimensional chain of Ising spins where exchange bonds occur with a probability 
decaying as the inverse power of the distance. Varying the power in this long-range model 
corresponds, in a one-to-one relationship, to change the dimension in spin-glass short-range 
models. Using different finite size scaling methods evidence for a spin-glass transition is found 
also for systems whose equivalent dimension is below the upper critical dimension at zero 
magnetic field. The application of a new method is discussed, that can be exported to systems 
in a magnetic field. 



1. Introduction 



Long-range (LR) models are such that their lower critical dimension is lower than 
the one of the corresponding short-range (SR) model. In particular, one can have a 
phase transition even in one dimensional systems, provided the range of interaction 
is large enough. One dimensional models with power-law decaying interactions 
actually allow to explore both LR and SR regimes by changing the power and 
enable to compare the ordered phase in and out of the range of validity of the 
mean field approximation. This is very useful for spin glass models that are known 
to have a rather complex states structure in the low-temperature phase in mean- 
field theory Whether this structure exists in finite-dimensional models with 
short-range interactions is, though, still a matter of debate. Theories alternative 
to the mean- field one have been proposed [1, 0], but SR systems are very tough to 
study analytically. Numerical simulations have been, thus, extensively employed, 
though with no conclusive indication on the nature of the spin-glass (SG) phase in 
finite dimension nor on the existence of a thermodynamic SG phase in presence of 
an external magnetic field. 

A one-dimensional spin-glass model with power-law decaying interactions 0- 
[SJl makes probes on larger linear sizes possible. Moreover, if bonds are diluted 
|lCH13l| the run time in numerical simulations grows simply as the size L of the 
system, rather than proportionally to L?' , as in fully connected systems. This is 
a fundamental issue because finite volume effects are strong in these models, so 
much that the very existence of the transition can be rather difficult to establish 
with canonical methods, e.g., when the interaction has a rapid decay and/or an 
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external magnetic field is applied. 



2. The ID Levy Spin-Glass Model 



The model investigated is a one dimensional chain of L Ising spins (uj 
Hamiltonian 



±1) and 



(1) 



The quenched random couplings Jij are independent and identically distributed 
random variables taking a non zero value with a probability decaying with the 
distance between spins cij and crj, rij = \i — j\ mod (L/2), as 



¥[Jij / 0] oc r 



for Tij 1 . 



(2) 



Non-zero couplings take value ±1 with equal probability. We use periodic bound- 
ary conditions and a z = 6 average coordination number. Links are generated by 
repeating zL/2 times the following process: choose randomly 2 spins at distance r 

with probability ?"^^/ X]j=i they are already connected, repeat the process, 

otherwise connect them. 

As the power p varies this model is known to display different statistical me- 
chanics behaviors [5r7]. For the diluted case ^JQ] they are reported in Tab. [U 
Equivalence between one-dimensional systems with power-law decaying interac- 
tions and D-dimensional systems with short-range (nearest neighbor) interactions 
can be approximately written as 12l |: 



1 



7/ 



D 



(3) 



where r] is the critical exponent of the space correlation function for the short-range 
model. The relationship is exact at the mean-field threshold p = 4/3 {D = 6, rj = 0) 
and approximated below. The analogue of a 3D spin-glass in zero magnetic field 
(with -0.384(9) < t] < -0.337(15) [li^il) would then be a system with p ~ 1.8. 
We will focus on this value in the present manuscript to introduce and test our 
method. Other model cases have been considered in Refs. id. Ill|. 



p 


D{p) 


transition type 


< 1 


oo 


Bethe lattice like 


]l:4/3] 


[6 : oo[ 


2nd Q^^gj,^ jYJp 


]4/3 : 2] 


[2.5 : 6[ 


2"^^ order, non-MF 


2 


2.5 


Kosterlitz-Thouless or T = 0-like 


> 2 


< 2.5 


none 



Table 1. From infinite range to short range behavior of the SG model defined in Eos Jll2t . 



Numerical simulations 



We simulate two replicas 0"^ 
The simulated sizes are L 



(1,2) 



using the parallel tempering algorithm (PT) [It 
2'^, with K, = 6,8,10,12. The interval between 
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3 

temperatures in the PT evolution is AT = 0.05. The number of samples is 
Nj = 6400 — 90000 depending on the size. All data used for our analysis are 
thermalized. Thermalization has been checked by measuring all observables on ex- 
ponentially growing time windows until the last two points coincide within the 
statistical error. 



4. Critical point 

The key observable to approach the critical behavior in disordered systems is the 
four-spins correlation function: 

c^4(x) = it(^FMMW5) (4) 
1=1 

and its Fourier transform C4{k). In order to determine the critical point, a corre- 
lation length-like observable is usually defined on the ID lattice as 0, [2o| 



2sin/ci/2 

with ki = 2-it/L. In Fig. [T]we present the ^/L curves whose crossing point should 
tend, as L — 7- oo, to Tc. In the inset we also show the behavior of 

XSG = LCM. (6) 

another finite size scaling (FSS) function for the present model (in which r] = 
'Hmf = p — ^ also for p > pmf)- Due to the statistical error it is not straightforward 
to identify clear crossing points for Moreover, in both above mentioned cases, 
to extrapolate a clear limit of Tc as L —)■ c» with a FSS interpolating function like 
a + bL~'^, cf. Fig. m three degrees of freedom are not enough and the interpolations 
are thus just indicative (see the following). 



^4(0) 



(5) 




0.6 0.7 



1.1 1.2 



Figure 1. FSS function ^/L vs. T for different simulated sizes at p = 1.8 for L = 26.8.10,12 i^seU XSG^''"^ 
vs. T for the same sizes. 
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Figure 2. Average (k) vs. [sin(fe/2)/7r]0-* for systems of size L = 2'5.8.10.12. Left: T = 0.9, shghtly 
below the critical region. Right: T = 1.3 > Tc- Insets: detail for low k (points) and comparison of the 
values in fe = with the fc — > limit of the interpolation A + _B[sin(fe/2)/7r]'''^ (lines). 



5. A novel method 

We can, otherwise, use the whole information provided by the C4^{k). In Fig. [2] we 
plot 1/C4,{k) vs. [sin(A;/2)/7r]^-i for p = 1.8 for ah simulated sizes both at temper- 
ature above Tc (right) and at a T ~ Tc (left). We observe that finite size effects act 
in opposite ways on the value of ^4(0) (i.e., xsg) and on the rest of the function 
6*4 (A: > 0) (see insets of Fig. [2|): while 1/(74(0) tends to its thermodynamic limit 
from above, 1/C4(A; > 0) and its interpolation at A; = tend to the thermody- 
namic limit from below. Even though 6*4(0) and lim^^o C'4(A;) are the same object 
for L ^ 00 their FSS scaling is qualitatively different. 
Interpolating 1/C4(A;) for small A; at a given size and temperature as 

F'^'ik) = AiL,T) + BiL,T)[smik/2)/TTf^ (7) 

we can analyze the L and T dependence of 

A{L, T) = lim C4{k) (8) 

k^O 

and determine the transition from the FSS analysis of the points at which 
A(L, T) = 0, rather than using the FSS of the crossing points of Eq. ([5]) (or 
C'4(0)L''~^). For p = 1.8, the behavior of A in L and T is plotted in Fig. [3j 
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Figure 3. Behavior of A{L,T), cf. Eq. l|8}, for L = 2^''^''^°'^^. The points at which ^ = are finite size 
estimates of Tc . 
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Figure 4. Comparison between the finite size estimates of Tc(L) by means of the canonical "crossings" 
methods (of the FSS functions g/L and XSG^''"^) and of the A{L,T) = method, cf. Eq. I©. The 
fit for Tc from A{L,T) = has a = 0.075 (full/red curve), while the other two interpolating curves 
(dashed/green for ^/L, dotted/blues for XSG^''"^) are only indicative (not enough d.o.f. to provide an 
estimate of statistical errors, no error bars defined on the Tc{L) from ^/L crossings). 



This method has the advantage of using high temperature data and one only 
needs to simulate systems down to the candidate Tc- As ^(oo, T) becomes negative, 
indeed, the functional form of the propagator in the paramagnetic phase breaks 
down and this provides evidence for a phase transition. In Fig. [J] we plot the finite 
size values of obtained by this method and we compare them with the estimates 
derived from FSS functions ^ and XSG- At /> = 1.8 we find Tc = 1.060(7) and 
l/u = 0.38(2) with a = 0.075 on the same data of the previous analysis (same 
statistics, same thermalization times, same sizes). 

We stress that also the definition of ^ as a correlation length in Eq. ^ is valid 
only in the paramagnetic phase and that below Tc this is just a scaling function 
without physical meaning. In that approach, though, in order to appreciate cross- 
ings of (T) I L curves at different L one has to simulate the system also at temper- 
atures below Tc, where thermalization times increase, and for at least five different 
sizes in order to provide enough points for a well defined FSS interpolation. 



6. Correlation length estimate 

We are interested in characterizing the above mentioned critical behavior by means 
of a growing correlation length, as it happens in ordinary continuous phase tran- 
sitions in systems without quenched disorder. Also for this analysis our starting 
point is the four-spins correlation function, this time in position space. As one can 
notice from Fig. [5l we can identify two different decays of the correlation at a given 
T > Tc, if we are able to study long enough Levy-glass chains. We observe a slower 
power-law decay as a; <C T, and a faster decay as x ~ L/2. Contrarily to 



what happens in short-range models at T* > 2.5 jl4l. l2ll|. this second decay is also 
power-law (with power equal to p) because the interaction correlation decays - by 
construction - as x~'^ and, therefore, the C4(x) cannot decay any faster. 

From the definition of r] as exponent of the anomalous decay of the correlation 
function as x~'^~^'^~'^, with r/ = 3 — p, we obtain the relation a + p = d+ l = 2 



22( 1 . We can, thus, think about interpolating the whole C4(x) behavior above the 
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Figure 5. Four spins correlation function at high temperature (T = 2 ~ 1.9Tc) in log-log scale for all 
simulated sizes L = 2®'*'!'''!^. Dotted lines: interpolation by means of Eq. Dashed-dotted line: x~>'. 
No finite size effects are present, apart from the values at i: ~ L/2, and a crossover from a power law with 
exponent a ~ 0.2 to a power law with p = 1.8 can be identified for L > 10. Inset: C4,{x) for L = 2^^ at 
T = 3 compared to x~^ *. 



1e+07 




1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 
T 

Figure 6. Correlation length parameter £ of interpolating function C4 * (x) vs. T in log scale for L = 2^ 
Inset: T behavior of fir parameters K and 5. 



critical point as 



(x) = Kx- 



1 + 



X\ 2<5(p- 
1) 



-1/5 



(9) 



with p = 1.8 and a = 2 — p = 0.2. 

From this interpolation we can look at the temperature behavior of the length- 
like parameter £, that is an estimate of the correlation length of the system, as far 
as the second power- law decay is observable (T > Tc). In Fig. [6] we display the 
behavior of ^(T) (and K and 6 in the inset) as T is lowered down to the estimated 
Tc- Fits with Eq. ([9]) are reasonable down to T = 1.1. 

As temperature approaches the critical value from above, the simulated systems 
are too small to appreciate the existence of a crossover length i ^ L and, thus, 
the C4(x) decay appears as a single power law. The same behavior remains below 
Tc, as shown in Fig. [71 This is incompatible with the onset of a plateau at any x 
whereas it is consistent with the clustering properties of the mean-field theory for 
spin-glasses. 
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Figure 7. Four spins correlation function across the critical region (T = 0.8,0.9,1,1.1) in log-log scale 
for all simulated sizes L = 2^'*'^'''^^ and relative interpolating functions. At all T, finite size effects are 
seemingly small. Dotted lines: interpolation by means of Eq. ((Oj at T = 1.1; for lower T the fit function 
is Ax~°' is used and no crossover to a power law with p = 1.8 is observed for any L, in agreement with 
£ ^ L. The T-dependence of the value of a is shown in the inset. 

7. Conclusions 

We have introduced a new FSS method to determine, by numerical simulations, 
the existence of a critical point in finite dimensional systems. Such method employs 
high temperature data. Thus, it requires lower thermalization times and disordered 
sample statistics with respect to canonical methods. It works well if one has a 
sufficient amount of points of the four-spins correlation function at low wavelength 
numbers k, that is, if the interaction range is not too broad (i.e., preferably out of 
mean-field) and linear size is long enough. 

We have tested the method in the case of a bond-diluted ID Levy spin-glass 
with a power-law decaying interaction outside the limit of validity of mean-field 
approximation and equivalent to a 3D nearest-neighbor interacting system on a 
cubic lattice. We have compared the results with those obtained by canonical FSS 
analysis on the same set of data and shown that they are compatible, cf. Fig. [H 

In position space, at T > Tc we have identified a crossover in the (power-law) 
relaxation decay of the four-spins correlation function from slow to fast (i.e., x~p) 
relaxation, cf. Eq. Q. The crossover takes place at a correlation length that grows 
exponentially as T decreases. At T < Tc the decay is well interpolated by a simple 
power-law, providing no evidence for the plateau predicted by droplet [1] and TNT 
theories \S. 
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